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Abstract. We discuss the formation and evolution of binaries which contain neutron 
stars or black holes. It is shown that in a stellar system which for 10 10 yr had star 
formation rate similar to the current one in the Galactic disc, the rate of neutron star 
mergers is ~ 2 x 10~ 5 yr _1 , consistent with the observational estimates. The rate of 
black hole and neutron star mergers is by an order of magnitude lower. The pairs 
of black holes form but do not merge because heavy mass loss by their progenitors 
efficiently moves the binary components apart. 

1 Introduction 

The wealth of observational features makes the merging of close neutron star 
(NS) binaries one of the most interesting phenomena in astrophysics. Apart 
from the expected bursts of gravitational waves (GWR, g), 7-ray bursts j| 
and production of r-process elements ]TT[ | are possibly related processes. Some 
NS binaries are observed as high mass binary pulsars (HMBP, Table |l|, see [^8| 
for references) , but these are possibly only a fraction of the entire population. 



Table 1. Observed population of high-mass binary pulsars 
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All estimates for the Galactic merger rates of NS binaries based on the 
observed population arrive at values smaller than ~ 10~ 5 yr -1 (e.g. Q, Q). 
Estimates based on simulations of the evolution of close binaries arrive at rates 
up to ~ 3 x lO^yr" 1 (e.g. Q). We will show that the estimated merger 
rates from observations and model computations may be reconciled. 
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Figure 1: Scenario of the evolution of massive close binaries. 



2 Formation of neutron star binaries 
2.1 Population synthesis 

"Theoretical" estimates for the birth and merger rates of NS binaries are usually 
derived from the models which combine the parameterized computations of the 
evolution of stars with prescriptions for the variation in the binary parameters 
due to the evolution of its components. By combining the model with distri- 
bution functions for the zero-age primary masses, mass-ratios, orbital periods 
and eccentricities, a complete picture of the population of binary stars is ob- 
tained. An example of the evolutionary scenario for a binary, which leads to the 



formation of two NS is given in Fig. W. 

Althou gh the concept of the formation of NS binaries is relatively well un- 
derstood (j^O), M), see however ||), the models imply some badly understood 
processes. The most crucial among these are the nonconservative mass exchange 
(especially the common envelope phase) and the supernova process. 
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A common envelope (CE) forms when the accretor in a binary is unable to 
acquire all matter lost by the donor. The outcome of the CE phase is either 
the coalescence of stars or the expulsion of the envelope. Models of the CE 
phase are still in an embryonic state. Therefore, a common envelope parameter 
a cc is usually introduced. The most common definition of a cc is the ratio of 
the binding energy of the CE to the orbital energy of the binary (e.g. Q). 
Models of the populations of close binary white dwarfs, low- and high-mass 
binary millisecond pulsars (see below) and some other objects may be fit to 
observations if a cc ~ 1 — 2. Though, formally, a cc ^ 1 means that sources other 
than orbital energy must be invoked in the expulsion of the CE, in fact, in such 
a parametric approach, this simply suggests that energy deposited into a CE 
has to be comparable to the orbital energy. Roughly, the ratio of semi-major 
axes of the orbit of a binary after and prior to the CE stage is af/ai tx a co , 
while the timescale for gravitational in-spiral is cx ah. This makes all merger 
rate estimates highly sensitive to a co . Estimates given below were found by 
Portegies Zwart & Yungelson |28) for a cc =2. 

A velocity kick may be imparted to the nascent NS. While the mechanism 
and distribution of natal kicks are still a matter of debate, there is firm ob- 
servational evidence for their occurrence Natal kicks may unbind binaries, 
which otherwise might have remain bound and, vice versa, conserve the binaries 
which without the kick would have been dissociated. An analysis of the role of 
kicks may be found e.g. in Q, Q . 

In |^8| we used the distribution for isotropic kick velocities v which has the 
functional form suggested in |23| ] with numerical parameters from Jl3[ : 

P(u)du — — ■ u = v/a, & = 600 km s _1 . (1) 

W 7T (1 + U 2 ) 2 ' W 

Implementation of Eq. (|l|) provides the best fit to the observed population of 
single pulsars close to the Sun (f|]). Peculiarity of this distribution is the 



combination of a peak at low velocity and of a high- velocity tail. Selection effects 
in the population of single radio pulsars, however, cause the exact shape of the 
kick velocity distribution to be ill constrained both at low and high velocities. 
Alternatives as e.g. a Maxwellian distribution with mean of 250 to 500 km s _1 
can be found in the recent literature (|l5[], p9|). 
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Figure 2: Probability distribution for orbital period and eccentricity for the 
population of NS binaries at the age T = 0, 10 s and 10 9 yr. The darkest shades 
correspond to the birthrate of 3.5 x 10 _7 yr^ 1 (top panel) and numbers of 110 
and 1100 (middle and bottom panels, respectively). Dots and asterisks show 
known HMBP. 

The results of computations without kicks, with a kick taken isotropically 
from Eq. [l], and with a kick from a Maxwellian distribution with a = 450 km s _1 
are presented in Table 2. Another parameters, like the initial mass-ratio distri- 
bution or initial orbital separation, affect the results by a factor of ~ 2. 



2.2 Merger rate of NS binaries in the Galactic disc 



A model which aims to predict a reliable birth and merger rate for NS binaries 
should provide also the characteristics of related types of binaries and HMBP. 
For the HMBP, the most useful observational parameters for comparison with 
the model are the orbital period P or b and eccentricity e. 

In |2£| the best fit between the observations of HMBP and the computations 
is obtained for a cc =2 and the distribution of kick velocities given by Eq. (|l]). 
Figure ^ shows probability distributions of the orbital parameters of the NS bi- 
naries with P or b < 10 4 day at an age of the population of T = 0, 10 8 and 10 9 yr. 
As demonstrated in [^8| , natal kicks "smear" P or b - e distribution. In the ab- 
sence of kicks only PSR B1913+16 appears in the region of the P or b - e diagram 
with high probability. Figure || suggests that the average age of the population 
of observed HMBP is several hundred Myr. Therefore it suggests that recycled 
pulsars are born with low magnetic fields. 

Table 3 compares the predictions of the preferred model from j2l| to the 
observations of NS binaries and related objects. The model is normalized to 
the current star formation rate in the Galactic disc of 4 M Q yr" 1 M min = 
0.1 M Q and 100% binarity are assumed. The adopted lower mass limit for the 
formation of a NS is 8M Q for an isolated star and 11.4 M Q for stars in close 
binaries. The model satisfactorily reproduces the occurrence rates of supernovae 
descending from massive stars, numbers of Be/X-ray and massive X-ray binaries 
(see also The model does not contradict observational estimates of the 

birthrate of single pulsars^] or the fraction of single recycled pulsars among all 
pulsars (see also [T^] ). It is important to notice that our computations do not 
violate the "Bailes criterion" Q , which limits the fraction of HMBP among the 
population of radio pulsars. 

The birthrate of NS binaries derived from computations is considerably 
higher than the observational one, which involves highly uncertain estimates of 
the incompleteness of the sample of observed objects |j| and their lifetimes fl33|] . 
Observational estimates are also likely to underestimate the actual birthrate, 
since a considerable fraction of old NSs in binaries may avoid recycling and 
after the death of the young pulsar never show up as pulsars[|. 

The results of computations agree with observations for the stages preceding 
formation of NS binaries, give a reasonable estimate of their birthrate and repro- 
duces their distribution in the P or b ~ e diagram. We therefore expect that com- 
putations also correctly predict the merger rate of NS binaries: 2.3 x 10~ 5 yr" 1 in 
the Galaxy. This estimate exceeds the most optimistic "observational" values 
by a factor of a few. However, about 10% to 20% of all NS binaries merge only 

6 These estimates strongly depends on the assumptions on the evolution of the magnetic 
field, the distribution of the interstellar material and the beaming fraction. For the latter 
/ = 0.3 was adopted. 

7 Lipunov et al. prf suggest that recycled pulsars constitute only several per cent of the 
total population of NS binaries; this estimate, however, weights heavily on the details of the 
accretion process on highly magnetized neutron stars. 



after the recycled pulsar has died in 0.5 - 1 Gyr (Fig. ||). Another ~ 10% merge 
within ~ 10 Myr, while the young pulsar is still visible. And we note that possi- 
bly only a fraction of old NS are recycled and show up in HMBP. This suggests 
a revision of the "observational" merger rate upward. 

Table 3. Descendants of massive binaries 
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Figure 3: Relative distribution of zero-age NS binaries over merger age. 



In |28| it is demonstrated that a higher average kick velocity in the supernova 
decreases the merger rate. In the most unfavored case of a Maxwellian kick 
combined with a mass ratio distribution peaked to and small a ce the merger 
rate may drop by an order of magnitude. 

The lifetime for a NS binary before it merges is typically below 1 to 2 Gyr 
(Fig. ||) . The merger rate is therefore dominated by the star formation rate in the 
last few Gyr. In massive (~ 10 11 M ) spiral (like the Milky Way) and irregular 
galaxies the star formation rate in the far past was probably significantly higher 
than at present |29|. However, past bursts of star formation contribute only 
little to the current merger rate of NS binaries. The same point has to be made 



with respect to elliptical galaxies, as most of their stars were formed in a short 
( ;$ 1 Gyr) burst; the overwhelming majority of NS binary mergers happened 
long ago. If we have two 10 Gyr old, equal mass galaxies with the same average 
star formation rate but in one case all stars are formed in a single initial burst 
and in the other the star formation was constant, the merger rate of NS binaries 
in the first galaxy is ~ 10 to 15% of that in the latter. 

Extrapolating our merger rates for NS binaries to the local Universe, using 
the blue luminosity of the Universe p3] and assuming h — 0.5, we obtain for the 
detection rate by the first generation of GWR detectors sensitive to distances of 
^ 25Mpc Q only ~ 1 in 200 yr! For advanced detectors (sensitive to a distance 
of 250 Mpc) the detection rate increases to ~ 5yr _1 . 

The merger rate for black hole + neutron star binaries (BH+NS) is expected 
to be about 20 times smaller than for NS binaries (28) and the rate of BH+BH 
mergers is virtually zero. The reason behind this is twofold. First, massive 
stars suffer heavy mass loss especially when they become Wolf-Rayet stars. 
E.g., evolutionary computations for solar composition stars with observationally 
inferred mass loss rates f3(| suggest that stars with an initial mass larger than ~ 
25 M Q have mass about 10 M Q prior to supernova. The stellar wind is expected 
to be isotropic and to carry away the specific angular momentum of the donor. 
Therefore mass loss in a stellar wind increases the orbital separation. Secondly, 
BH are not expected to receive sizable kicks which may, if "properly" directed, 
decrease the binary separation. 

The signal-to-noise ratio for gravitational waves detectors is given by B : 

S/N oc (M 1 M 2 ) 1 /2(M 1 + M 2 )-^ e D-\ (2) 

where M 12 are the masses of coalescing objects and D is the distance to them. 
For M BH =4 4- 17 M suggested by observed BH candidates, one would expect 
that a BH+NS merger is detectable to about 1.5 to 2.5 times as far as NS+NS 
merger. Hence, the rates of detections of NS+NS and BH+NS mergers may be 
comparable. 



2.3 Merger rate of NS binaries in globular clusters 

Also in the dense cores of globular clusters (GC) mergers between NS may occur. 
The fraction of primordial binaries which evolve to the stage where GWR drives 
the components into coalescence is, however, small. The majority of them are 
ejected from the cluster or dissociated by the encounter with another star. Some 
of the ejected binaries still result in a NS binary and finally merge. Single NS in 
the cores of high-density GC may find another star and be captured by it either 
in a two- or three body encounter. If two NS are able to form a binary in such 
a way this may increase the NS binary merger rate for post-collapse clusters 
considerably. 

A reliable estimate for the occurrence rate of such mergers is not trivial as 
the evolution of the binaries and the stars are coupled in a complex fashion. 



For solving this problem direct N-body simulations with a realistic number of 
stars are required, for which the hardware simply does not exist at present. Re- 
cently, however, Portegies Zwart et al. ]26| ], |27j developed a simplified model 
for simulating the interaction between evolving stars and binaries while taking 
dynamical encounters into account. The microscopic dynamics of close encoun- 
ters between single stars and binaries are solved using semi-analytic and 3-body 
techniques, while the macroscopic evolution of the stellar system is assumed to 
be unaffected by its slowly changing stellar population. 

The results of these computations concerning the merging of NS binaries 
may be summarized as follows: In the model for a core-collapsed cluster, the 
fraction of merging primordial binaries is ~ 60%. The fraction of NS binaries 
among these is ~ 0.004. The time interval between core collapse and Hubble 
time is ~ 5 Gyr. If one takes a globular star cluster with 10 6 stars of which 10% 
are in the region with a high density and 25% primordial binaries, one obtains 
a rate for NS binary mergers of ~ 1.2 x 10 -8 yr -1 . 

Of the binaries which are formed by tidal capture, which occurs once every 
10 7 years, the merger rate for NS binaries is only 4 x 10~ 10 yr _1 (almost all 
tidally formed binaries coalesce, but only a fraction of 0.004 contains two NS). 

For a total of ~ 250 GC in our Galaxy merger rate for NS binaries becomes 
~ 10 _7 yr _:L . The detection rate for the early GWR detectors is therefore only 
^ 1/1000 years (depending on the fraction of dense star clusters). Thus, the 
rate of NS mergers in the field under any circumstances is much higher than the 
rate of mergers in GC. 

3 Conclusion 

"Theoretical" and "observational" merger rates of NS binaries produce similar 
results if the ejection of CE is highly efficient and nascent neutron stars receive, 
on average, velocity kicks of the order of (200 -i- 300) km s _1 . 

The rate of NS mergers is ~ 2 x 10 -5 yr _1 for a galaxy with constant astration 
rate of 4 M Q yr _1 (and M min = 0.1 M Q and 100% binarity). This rate is 
uncertain by a factor of a few mainly due to uncertainties in the kicks velocity 
distribution and the efficiency of energy deposition into the CE. 

Merger rate is mainly determined by the current star formation rate and 
much less by the star formation history of the Galaxy. The time to coalesce two 
NS stars is typically a few times 10 8 years. 

Extrapolating the Galactic merger rate of NS binaries, we arrive at a detec- 
tion rate of once every ~ 200 yr for the first generation GWR detectors (which 
are expected to be sensitive up to ~ 25Mpc). BH+NS mergers may be regis- 
tered at the rate comparable with binary NS mergers. 

Mergers of BH+BH binaries are not to be expected because the severe mass 
loss in the stellar winds of their very massive progenitors results in too large for 
merger orbital separations. 



Only ~ 1% of the GWR bursts are expected to originate from GC. 
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